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The  Defense  Communications  Engineering  Center  (DCEC)  Technical  Notes 
(TN's)are  published  to  inform  interested  members  of  the  defense  community 
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I.  INTRODUCTION 


A military  telecommunications  technical  control  facility  is  composed  of 


a group  of  technical  control  operators  who  have  the  responsibility  of  re- 
storing circuits.  When  one  of  the  circuits  goes  down  a technical  controller 
is  assigned  to  it  to  ensure  it  is  repaired  and  brought  up  again.  In  a recent 
study  [1],  it  was  pointed  out  that  a technical  control  facility  can  be 
considered  as  a multiserver  queueing  system.  Furthermore,  one  possible  op- 
erating rule  for  the  facility  had  not  previously  been  considered  in  the  queue- 
ing theory  literature.  It  is  the  purpose  of  this  paper  to  present  a queueing 
model  for  this  rule. 

The  rule  can  be  stated  as  follows:  Independent  of  the  number  of  circuits 

requiring  corrective  action,  one  technical  controller  periodically  leaves  the 
facility  to  perform  other  work.  After  he  is  finished  he  returns  to  the  facil- 
ity. Although  the  analysis  of  this  paper  assumes  onlycne  of  the  technical  con- 
trollers may  leave,  the  extension  to  the  situations  when  more  than  one  may  leave 
follows  via  the  same  sort  of  analysis. 

In  queueing  theory  terminology,  the  customers  are  the  circuits  and  the 
servers  are  the  technical  controllers  who  are  assigned  to  the  facility.  For 
the  particular  operating  rule  considered  above,  we  will  be  considering  an  s 
server  system  where  periodically  one  server  leaves  the  facility  to  return 
sometime  in  the  future.  In  queueing  theory  literature,  we  are  considering  a 
system  where  one  of  the  servers  may  break  down. 

The  analysis  of  queueing  systems  with  breakdowns  has  been  considered  in  the 
past.  Gaver  [2]  introduced  and  used  the  notion  of  completion  times  to  analyze 
the  case  where  the  number  of  servers  equals  one.  Several  other  related  papers 


* 


have  appeared,  e.g.  Jaiswal  [3],  We  will  not  mention  them  all,  although  one,  re- 
ference [4],  should  be  noted  as  it  was  the  first  paper  to  consider  the 
multiserver  case.  In  that  paper,  Metrany  and  Avi-Itzhak  considered  a system 
where  all  the  servers  may  go  away  or  break  down  The  operating  rule  under 
consideration  in  this  paper  is  sort  of  a hybrid  of  the  rules  considered  ^'n 
[2]  and  [4].  That  is,  only  one  of  the  s servers  leaves  the  facility  or  breaks 
down. 

In  this  paper,  we  give  an  exact  analysis  of  the  case  where  the  circuits 
(customers)  are  assumed  to  break  down  in  accordance  with  a Poisson  process 
and  the  length  of  time  required  to  fix  the  circuit  has  an  exponential  distribu- 
tion. There  are  s ( >1 ) technical  controllers  (servers)  in  the  facility  and 
one  of  them  periodically  leaves  the  facility.  We  assume  that  the  length  of  time  he 
leaves  and  stays,  each  has  an  exponential  distribution.  In  section  II  we  give 
an  exact  analysis  of  this  queueing  system.  Some  numerical  examples  are  given 
in  section  III, as  well  as  several  possible  approximations.  Finally,  section 
IV  contains  some  concluding  remarks. 
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II.  ANALYSIS 


In  this  section,  we  give  a mathematical  analysis  of  the  queueing  system 
described  in  section  I.  Thece  are  s servers  and  periodically  one  of  the 
servers  leaves  the  system  to  return  at  some  future- time.  We  assume  that 
customers  arrive  to  the  system  in  accordance  with  a Poisson  process  with' 
parameter  x.  The  service  time  of  an  arriving  customer  has  an  exponential 
distribution  with  mean  u""*.  For  the  server  who  periodically  leaves  the  sys- 
tem, we  assume  the  random  variables  representing  length  of  time  he  remains  in 
or  out  of  the  system  have  an  exponential  distribution  with  means  and 
respectively.  All  random  variables  are  assumed  to  be  mutually  independent 
and  we  assume  the  system  has  an  infinite  waiting  room. 

Let,  for  n=0,l,2,...,  and  i=0,l 

pn,i  = Pr(Q=n,  Y=i},  (1) 

where  Q is  the  steady  state  number  of  customers  in  the  system  and 

0 if  the  server  away  from  system 

1 if  the  server  in  system. 

If  a customer  is  receiving  service  when  the  server  decides  to  leave,  the  cus- 
tomer is  returned  to  the  head  of  the  queue-  Waiting  customers  are  serviced 
on  a first-come,  first-served  basis. 

The  steady  state  equations  for  P . are 
...  n , l 

(*+nu+e)Pnp  * xPn-i,o+(n+1)vjPn+l,0+aPn,l  ; 

(A+(s-l)p+e)P^0  = *pn_1>o+(s-l)ppn-H,0+0‘Pn,l  ; 


n<s-2 


n>s-l 


(2) 
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(x+nu+a)PnJ  = APn_ltl+(n+l),iPn+1J+3Pn>0  ; "<s-l 
(X+Sp+a)Pn>i  = APn_1J+SpPn+1J+BPnj0  ; n>S  . 

where  P_-|  g = -1  1 ~ Let,  f°r  an<^  Mil* 

PjU)  - Z p„  ,z". 

1 n=0  n>1 

Multiplying  equations  (2)  and  (3)  by  zn  and  combining  one  gets 


(-xz*+(B+x+(s-l)u)z-(s-l )n)PQ(z)-azP1 (z) 


=u(l-z)  Z (n+1 -s)P  nzn 
n=0  n,u 


7 3>” 1 

(-XZ  +(a+X+Sp)z-Sy)P, (z)-6zPn(z)  = u(l-z)Z  (n-s)Pn  ,Z. 

' u n=0 

These  equations  can  be  written  in  the  following  matrix  equation: 


A(z)P(z)  - B(z) 


where 


P(z)  = 


P0(Z) 

P,(*> 


B(z)  = 


w(l-z)  Z (n+1 -s  )P  nzn 
n=0  n,u 


u(l-z)  Z (n-s)P  ,zr 
n=0  n’’ 


an(z)  -oz 
A(z)  = 0 

-6z  a^z) 

• 2 2 
with  ag(z)=-Xz  +(e+X+(s-l )y )z- (s-1 )y  and  a^(z)=-Xz  +(a+x+sy)z-su. 
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If  det(A(z))  is  the  determinant  of  A(z)  then  the  solution  for  P^(z), 
1=0,1,  is  given  by 


s-1 

X 


s-1 

n.  „ X 


P0(z) 


^1-zJ^l<z)n=0<n+1-s)Pn,0zn+azn=0(n-s>Pn,l2n}  (8) 


det(A(z) ) 


and 


s-1  s-1 

P (2)  = t,(1-z>(a0U)n^0<n-5)|’n,lA8znS)<nt1-s)Pn,0zni  ■ (9) 

' det(A(z) ) 


where 


det(A(z) ) = aQ(z)a1 (z)-aBz2  . 


(10) 


Thus,  the  solution  depends  on  P . (i=0,l;  n=0,l  ,2, . . . ,s-2)  and  P , , 

n $ i s — l j I • 

These  probabilities  may  be  expressed  in  terms  of  Pg  g and  Pq-j.  Consider  equa- 
tions (2)  and  (3),  for  n=0,l ,2, . . . ,s-2  we  have 


Pn,0  = C0(n)p00+Cl(n)p 


01 


and  for  n=0,l ,2, . . . ,s-l 


where 


and 


Pn,l  D0^P00+[M">P01 

nuCg(n)  = (x+6+(n-l)u)Cg(n-l )-xCg(n-2)-aDg(n-l ) 
npC i ( n ) = (A+B+(n-l  )y  )C-|  (n-1  )-xC-|  (n-2)-aD-|  (n-1 ) 

npDg(n)  = (x+a+(n-l  )u)Dg(r,-l  )-A0g(n-2)-BCg(n-l ) 
npD-|(n)  = (x+a+(n-l  )m)D-|  (n-1  )-xD-|  (n-2)-3C^  (n-1 ) 


(ID 

(12) 


(13) 


(14) 
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with  Cq(0)=D-|  (0)=1 , C-|  (0)=Dg(0)=0  and  for  i=0,1  (n)  = D^  (n)=0  when  n<0.  This 

shows  that  the  solution  to  the  problem  rests  on  determining  the  two  unknowns 


P00  and  P0T 


One  equation  in  these  unknowns  can  be  found  by  the  normalizing  condition*, 

co  i 


that  is.  2 2 P • = 1 . From  equations  (8)  and  (9),  using  L'Hospital1^  rule, 

n=0  i=0  11,1 


one  gets  (p=x/p) 


s-1 


s-p-a/(a+g)  = n50^C(s-1-n)Cg(n)+(s-n)D0(n)]P00 


(15) 


+ [(s-l-n)C1(n)+(s-n)D1(n)]P01}. 


We  need  to  generate  another  independent  equation  in  Pqq  and  Pg-| 

Let  us  consider  the  equation  det(A(z))=0  by  using  Sturm  sequences,  [5]; 
one  can  show  there  are  exactly  two  roots  of  this  equation  in  [0,1].  Obviously, 
one  of  the  roots  is  equal  to  one,  and  note  that  if  s=l , the  other  is  equal  to 
zero.  Thus,  there  is  exactly  one  root  of  det(A(z))=0  in  [0,1);  we  denote  it 
by  Zg.  Returning  to  equations  (8)  or  (g);  the  numerator  of  each  of  these 
equations  must  vanish  at  Zg.  This  gives  us  two  more  equations  in  Pgp  and  Pg-j . 

A little  algebra  will  show  that  they  are  equivalent,  and  so  we  have 


s-1  n 

0 = Pgp{  ^[a-| (zQ)(n+l-s)Cg(n)-aZg(s-n)Dg(n)]Zg} 

s-1  n 

+ PQ1{  I^[(a1(Zg)(n+l-s)C1(n)-aZQ(s-n)D1(n)]Zg} 


06) 


From  equations  (15)  and  (16)  we  can  now  give  an  expression  for  Pqq 


and  P 


qr 


s-1 

2 [(s-l-n)C,(n)+(s-n)D.(n)] 


For  i=0,l , let 


(17) 


d s-1  n 

ni  = I^[a1 (z0)(n+l-s)Ci(n)-az0(s-n)Di(n)]z0  (18) 

then 

- (s-p-ot/ (a+3) )n-| 

P0.0  = 

YOnl ~Y1 n0 

% 

and 

(s-p-a/(a+g) )nn 

Pm  = • (20) 

Although  no  simple  expression  exists  for  the  expected  number  of  customers 
in  the  system,  E(Q),  one  may  differentiate  equations  (8)  and  (9)  to  obtain 
the  desired  result.  For  completeness,  we  include  the  resulting  expression 


E(Q)  = (u/0T ){ [©n(cs+&-X+(s-l)y)-0, (a+e)]  2 (n-s)P  , 
i u n=0  n’ 


+ [0n(a+g-X+su)-e,(a+e)]  I (n+1 -s)P  n+0n(a+3)  - n(n-s)Pn  , 
u 1 n=0  n,u  u n=0  0,1 


+ 0n(a+e)  £ n(n+l-s)P  n> 
u n=0  n,u 


where 


0q  = B(A-Sp)+a( A- ( s-1 )p) 

0-|  = (X— 3— ( S-1 )w) (a-X+Sp)+BX+aX+af, 
and  the  P .'s  are  found  from  (11),  (12),  (19)  and  (20). 

n , i 
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III.  SOME  NUMERICAL  EXAMPLES  AND  APPROXIMATIONS 


In  this  section  we  give  some  numerical  examples  of  the  results  found 
in  section  II,  as  well  as  discussing  several  possible  approximations  to 
the  system.  Table  I gives  a comparison  of  the  system  for  the  same  traffic 
intensity,  p=X/p(y=l),  and  two  different  values  of  s,  (s=5  and  s=l 0) . We  note 
that  as  g-*0  the  system  behaves  as  ans-1  server  system. 

TABLE  I.  EXPECTED  NUMBER  OF  CUSTOMERS  IN  SYSTEM  (E(Q)) 


From  this  table  one  can  see  that  the  system  with  5 servers  is  more 
sensitive  to  one  of  the  servers  leaving  than  the  system  with  10  servers. 


i 


Of  course,  this  was  to  be  expected  since  when  one  server  leaves  in  the  five 
server  system,  the  total  system  capacity  is  reduced  by  20 %,  whereas,  in  the 
other  case  only  by  10%. 

Although  a straightforward,  standard  analysis  was  presented  in  section  II 
development  of  some  simple  approximations  to  the  system  that  do 
not  require  finding  Zq  may  be  desirable.  One  approach  may  be  to  tend  the  de- 
sired measure  of  performance  (e.g.,  expected  queue  length,  average  waiting  time) 
for  an  s and  s-1  server  system,  and  then  take  the  convex  combination  of  these 
measures  of  performance  based  on  the  proportion  of  time  the  system  is  operating 
like  an  s server  (e/(a+&))  and  an  s-1  server  system  (ct/(a+8)).  This  approach 
is  not  very  good,  since  the  traffic  intensity  may  be  such  that  for  an  s-1 
server  system,  the  system  blows  up.  For  example,  suppose  s=5,  a=l,  0=2,  p=l, 
and  A=4;  since  p=4,  the  measure  of  performance,  say  expected  number  in  the 
system,  for  a 4 server  system  would  be  infinite.  Thus,  taking  a convex  com- 
bination of  this  number  and  the  case  for  s=5  would  result  in  an  infinite  number 
of  customers  in  the  system.  Of  course,  this  is  wrong  since  p<s-a/(a+8),  and 
thus  the  Cesaro  mean  converges. 

A more  promising  approximation  is  to  treat  the  system  as  a non-integer 
number  of  servers.  That  is,  consider  the  system  to  be  an  M/M/s'  system  when 
s'=s-a/(a+8).  The  only  problem  with  this  approach  is  in  computing  the  desired 
measure  of  performance  for  a non-integer  number  of  servers.  From  [7]  one  can 


express  most  of  the  measures  of  performance  in  terms  of  the  Erlang  Loss 
Formula,  E(s,p),  where 


Thus,  the  only  problem  is  to  determine  Erlang's  Loss  Formula  for  a non-integer 
number  of  servers.  As  suggested  in  [8]  we  use  the  following  extrapolation 
formula  for  this  quantity, 

E(s',p)  = E([s'],p)[-rslj+1+pE([svjjp)j  j (22) 

when  [x]  is  the  greatest  integer  less  than  or  equal  to  x.  Tables  II  and  III 
give  a comparison  of  this  approximation  with  the  exact  results.  The  results 
are  presented  for  s=5  and  s=10,  and  various  values  of  a and  6.  We  note  that 
for  the  case  where  a=0  or  6=0  the  approximation  is  exact. 


TABLE  II.  COMPARISON  OF  APPROXIMATION  FOR  E(Q)  WITH  s=5 
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Two  points  are  immediately  discernible  from  these  tables.  First, the  approximation 


F 


always  underestimates  the  exact  results,  and  second,  the  approximation  is 
better  for  the  larger  number  of  servers.  In  both  cases,  the  approximation 
is  extremely  good  and  seems  to  be  worst,  in  a relative  sense,  when  a=$. 

Although  this  approximation  is  extremely  good,  the  only  real  computational 
problem  in  the  exact  analysis  is  in  determining  Zq, the  root  of  det(A(z))=0 
inside  of  [0,1].  A numerical  investigation  of  the  relationship  between  PqQ 
and  Pqi  revealed  a very  interesting  fact;  that  is, 

ePqo  * aP0>l 

where  - means  approximately.  Table  IV  gives  some  typical  results  that  were 
used  in  generating  this  observation 

TABLE  IV.  NUMERICAL  INVESTIGATION  OF  THE  EQUATION  SP00=cxP01  (s=5) 


p 

P 1/3  P 

*00  K0.1 

pao  1 P0,l 

P 3/^  P 

kQ0  kQ1 

1 

.091928  .27582 

.18379  .18380 

.27561  .09187 

Vy 

.3333 

.9999 

2.9999 

2 

.03338  .10017 

.06626  .06636 

.09872  .03292 

Vo; 

.3333 

.9986 

2.9992 

3 

.01125  .03380 

.02146  .02160 

.03056  .01020 

Vo; 

.3329 

.9933 

2.9957 

3.9 

.00326  .00981 

.00518  .00527 

.00534  .0017877 

.3322 

.9825 

2.9878 

This  relationship  suggests  another  possible  approximation;  that  is,  to 
use  equations  (15)  and  (23)  to  solve  for  PQp  and  directly  without  having 
to  find  Zq.  Once  Pqq  and  Pg-^are  found,  Pn  1.(n=l,2,...  and  i=0,l)  can  be 
determined  recurs ively,  or  expected  value  measures  of  performance  can  be  ob- 
tained as  was  done  in  the  case  of  E(Q)  in  section  II.  We  note  that  this 
approximation  has  two  advantages  over  the  other  one:  it  is  more  accurate,  and  it 
can  be  used  to  generate  the  complete  probability  distribution  if  required. 

The  greater  accuracy  is  shown  in  Table  V.  In  that  table,  approximation 

I is  the  one  using  the  non-integer  number  of  servers  and  approximation  II,  the 
one  just  suggested. 

TABLE  V.  COMPARISON  OF  TWO  APPROXIMATIONS  FOR  E(Q)  (s=5) 


^^^oi/(a+B) 

p/S 

.25 

.5 

cn 

APPROX.  I 
.2  APPROX.  II 
EXACT 

1 .0015 
1.0021 
1.0022 

1.0025 

1.0036 

1.0037 

.3 

[|EE| 

m 

1.5292 
1.5335 
1 . 5338 

.4 

■ : ■ 

2.0813 

2.0956 

2.0965 

HI 

.5 

2.6811 

2.6973 

2.6987 

2.7549 

2.7860 

2.7881 

2.8649 

2.8890 

2.8910 

.6 

3.4907 

3.5198 

3.5222 

3.6944 

3.7540 

3.7577 

4.0104 

4.0504 

4.0548 

.7 

m 

5.3544 

5.4664 

5.4721 

6.4231 

6.5050 

6.5106 
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IV.  CONCLUSIONS 

I 

In  this  paper  we  have  presented  a queueing  model  of  one  possible  operating 

rule  of  a techncial  control  facility.  The  numerical  examples  considered  in 

section  III  pointed  out  serveral  interesting  facts.  First,  the  effect  of  one 

of  the  servers  leaving  the  facility  for  periods  of  time  is  more  pronounced 
' 

when  the  number  of  servers  is  small.  Second,  the  two  unknowns,  Pg  g and 
Pg  are  approximately  related  as  gPQ  g=PQ  ^a.  This  fact  can  be  used  to  give 
one  possible  approximation  to  the  system.  Another  possible  approximation, 
investigated  in  section  III,  is  to  consider  an  equivalent  system  when  there 
is  a non-integer  number  of  servers,  s'  with  s'=  -a/(a+g). 

It  is  hoped  that  the  exact,  as  well  as  approximate,  analysis  of  a possible 
operating  rule  for  a technical  control  facility  may  provide  some  insights  into 
the  operation  of  such  a facility  as  well  as  being  used  in  performing  some 
trade-off  studies  for  the  system. 
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